#ifndef NUTATION_H
#define NUTATION_H

#include <math.h>
#include <wchar.h>  // 宽字符支持

#define rad (180 * 3600 / M_PI)  // 每弧度的角秒数

// IAU2000B 章动序列
extern const double nuTab[];

// 中精度章动计算表
extern const double nutB[];

// 章动计算
double *nutation(double t, double zq);

// 计算赤经章动及赤纬章动
double *CDnutation(const double z[3], double E, double dL, double dE);

// 中精度章动计算
double *nutation2(double t);

// 只计算黄经章动
double nutationLon2(double t);

#endif // NUTATION_H